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ABSTRACT 

This document describes the computer program developed to describe the 
three-dimensional motion of the Dynamic Docking Test System (DDTS) active 
table. The input consists of inertia and geometry data, actuator struc- 
tural data, forcing function data, hydraulics data, servo electronics 
data, and integration control data. The output consists of table responses, 
actuator bending responses, and actuator responses. 
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1.0 GENERAL INFORMATION 

1.1 COMPUTER SYSTEM 

This program was written in FORTRAN V language for use on the UNIVAC 1108 
computer with the EXEC II operating system. The program can easily be 
converted to the EXEC 8 operating system. The plotted output is done on 
the Stromberg Datagraphix's SD-4060 microfilm plotter. 

r 

1.2 PURPOSE 

The program simulates the motions of the NASA JSC Dynamic Docking Test 
System (DDTS) active table. Given a description of the table mass and 
geometry, the actuators, the hydraulic system, the electronics, and the 
forcing function, the program outputs time-histories of table responses, 
actuator bending data, and actuator responses. Responses to the following 
input forcing functions are calculated: ' _ 

a. Step velocity conmand 

b. Sinusoidal position commands 

c. Step external fore*-.- on table c.g. 

d. Sinusoidal external force on table,. c.g. 

1.3 LIMITATIONS 

Forcing functions are limited to those listed in Paragraph 1.2. Dynamics 
of the simulator due to docking are not modeled. Actuator control system 
components are limited to those shown in Figure 4. In the actual DOTS, 
there are notch filters in both the velocity command line and the forward 
loop which are not included in this simulation. 

i. 

2.0 PROCEDURE ' 

2.1 PROGRAM NAME 3 

The program acronym is NADS from NASA Advanced Docking ^stem. 


'-s' 
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2.2 NOMENCLATURE 


Nomenclature for NADS - 


Program 

Engineering 

Description 

A 

tA] 

Transformation matrix from table to inertial 
coordinates 

A 

Aj, Ag 

“Push" and "pull" stroke working areas of 
actuators 

AA 

•« 

’p 

Actuator acceleration 

AL 

V 

Actuator length 

ALPHA 

a 

Break frequency of first order filter 

AL3D 

••• 

Ip 

Actuator jerk 

AV 

’p 

Actuator velocity 

BETA 

3 

Break frequency of first order filter 

BETAE 

^e 

Equivalent hydraulic system bulk modulus 

BP 


Viscous damping coefficient of actuator 

C 


Column of generalized forces for equations of 
motion solution 

CAA 

• • 

’c 

Commanded actuator acceleration 

CAL 

’c 

Commanded actuator length 

CAV 

# 

’c 

Conmanded actuator velocity 

CP 

s 

Leakage coefficient across piston seals 

DLTAI 

AX^, &yj, 
AZj, A0. 

Ai|)» A(J) 

Sinusoidal amplitudes of translational commands 
for table c.g. and of table Euler angles 

EIR 

“r 

Bending modulus of piston rod 

FF 


Coulomb friction force of actuator 

FH 

fH 

Total hydraulic and friction forces acting 
on pistons ^ 



f 


D2-1185M-2 


2.2 (Continued) 


Program 

FMEXT 

FP 

FRQNCY 

lAC 

IFIRST 

IM 

INDKTR 

INER 


IPLOPT 

IPROPT 

IXF 

KC 

KF 

KG 

KPF 

KR 

KRC 

LC 


Nomenclature for NADS (continued) 
Engineering Description 

External forces and moments 
Net forces on actuator piston 

r 

u) Command signal frequency 


'"ext* ^EXT 

Fn 


I 


AC 


I 


I 


XX ’ W’ 

zz’ ^xy’ 
'yz 


XZ’ 




rc 


c 


Mass moment of inertia of cylinder (excluding 
the mass of the piston) about floor swivel 
joint 

Initialization indicator for mass matrix and 
geometry 

Mass matrix and geometry update option indicator 

Stroking or matrix inversion error indicator 
to terminate the integration process 

Moments and products of inertia 


Plot option indicator 

Print option indicator 

External force and moment option indicator 

Valve pressure flow coefficient 

Displacement feedback and command gain 

Electronics and valve forward loop gain 

Pressure feedback loop gain 

Velodty feedback loop gain 

Velocity command gain 

Distance from floor swivel to center line of 
piston rod seal at end pf cy]inder 
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2.2 (Continued) 


Nomenclature for NADS v continued) 


Program Engineering 


Description 


Maximum stroke of actuators 


Length of piston rod 

Retracted length (between swivel joints) of 
actuators 


Mass matrix and mass matrix inverse 

Moment acting about table c.g. from hydraulic 
and friction forces 


Effective rigid lateral mass of actuator 
assembly 

Mass of piston rod and piston 

Effective bending mass lumped at rod seal of 
cylinder 

Table mass 


NFREQ 

NFFREQ 


Number of table displacement frequency cases 

Number of external force and moment frequency 
cases / 


NPLTS 

OMEGA 

OMEGAC 

OMEGAE 


OMEGAF 


OMEGAS 


OMEGAV 

OMEGPF 


(Uj, 


Vl’ V2 


Number of plotted time points 

Break frequencies of first order filters 

Displacement command signal frequency 

Actuator bending frequency 

Frequency of sinusoidal external forces and 
moments • 

Frequency of second order filter on displacement 
and v^l oci ty feedbacks 

Frequency of valve dynamics 

Break frequencies of pressure feedback filters 
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2.2 (Continued) 


Nomenclature for NADS (continued) 


Program 

Engineering 

Description 

OUTFRQ 

- 


Output frequency for printing and plotting 

PS 



Supply pressure 

RS 

"s 


Inertial vector components of actuator length 

RXA 

"xa 


X axis table station of actuator swivel joints 
with respect to the table c.g. 

RYZA 

••ya- 

■■za 

Y, Z table coordinates of swivel joints with 
respect to the table c.g. 

T 

[T] 


Transformation matrix transforming vectors from 
table coordinates to local actuator coordinates 

TC6C0 



Initial inertial coordinates of table c.g. 





TDIC- 

r-c: 

{M 

Time dependent inertial commands 

> 


1%; 

M“-c) 

► 

TEAO 

®0’ '*'0’ *^0 

Initial Euler angles of the table coordinate 
system with respect to the inertial system 

TEND 

- 


Last integration time 

TIME 

t 


Time ' 

TITLE 



Title to be printed at top of first page of 
output 

TPLOT 



Time point at. or after which output for plots 
is made ' 

TPRINT 

- 


Time point at or after which printed output 
Is made 

/ , 

TSTART 


- 

First integration time v 
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2.2 (Continued) 


Program 
VO 

X 

XDOT 
XO 
YZF 

ZETAE 
ZETAS 

ZETAV 

2.3 METHOD 

This section contains a brief description of the physical system for which 
the program was written and the mathematical equations used to describe 
the motion of the system. The equations are described in detail in the 
referenced document. 

The DOTS active table is a triangular platform supported by six hydraulic 
actuators as shown in Figure 1. The table is capable of six-degree-of- 
freedom motion controlled by the six actuators. 

Three coordinate systems are ’used to describe the motion of the table and 
actuators. These coordinate systems are shown in, Figure 2. The inertial 
coordinate system origin is on the simulator centerline in the plane of the 
floor swivel joints. Table motion commands and responses are expressed in 
the inertial coordinate system. Table coordinates are body fixed coordi- 
nates whose origin is at the table center of gravity. Actuator coordinates 
are used to describe actuator motions. 


Nomenclature for NADS (concluded) 

Engineering Description 

V Initial hydraulic volumes of push and pull 

strokes of fully retracted actuator 

Variable array (output by integration procedure) 

Derivative array 

Xq Initial condition array 

y^, Zx Y and Z inertial coordinates of floor swivel 

joints 

Cg Damping constant for actuator bending 

r Damping constant of second order filter on 

displacement and velocity feedbacks 

C Damping constant of valve dynamics 






Figure 1. DOTS SlmuUto.' Facility 
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TABLE 

COORDINATES 



Xj ACTUATOR 

COORDINATES 



Figure 2. Active Table Coordinate Systems 
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2.3 (Continued) 

Table equations of motion are written in the body fixed table coordinates 
as follows: 


/ • 



IxI = [m]*Mc| (1) 

where: {xl is a‘ column of* accelerations for each degree of freedom (six 
degrees of freedom for the table and two elastic degrees of 
freedom for each actuator) 

[M] is the 18 x 18 coupled mass matrix 

|C| is a column of general i 2 ed forces for each degree of freedom 

The mass coupling effects of the actuators due to table motions are 
determined by Lagrange's method. The three-dimensional rigid motions of 
the actuators are completely constrained (i.e. , they are dependent upon 
the table motions). The mass matrix is shown in Figure 3 in upper tri- 
angular form. 


The, column of generalized forces includes the velocity terms in the 
equations of motion, the total forces exerted on the table by the actua- 
tors, externally applied forces and moments, and actuator bending stiffness 
and damping. 
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2.3 (Cpnijnued) 

The Euler angles e» \|», are used to transform velocities in the table 
coordinate system to the inertial coordinate system 


I - J : 


i ' , ; ’ \ 


6 

I 


\ I 




1 

^1 

= [A] I 


• 

1 

) 


. sioiii 
COS\j; 


■q cosa 

COSil) 

0 sin4» cos4« 

1 -coS({itan4> sin^itan^j) 




b) 


(t) 

V z 


(3) 


(4) 


v/here : 


1 [A] " 


C6*C\j) -C4i*C6*Sii;+S0*Sij> S(|)*C6»Si}i+C<i>*Se 


Sij) 


C({.»Ct!) 


-S(j)*C4» 


C$*S6*S\|»+S(j)*C0 “S({i*S0 *Si}i+C4)*C0 J 


(5) 


! C = cosine 

t 

S = sine 

' ' I 

Ea-ch actuator is modeled as a flexible rod with pinned ends and is free 
to bend in its firs, lateral mode in two orthogonal directions. 

,,, , 1 

Hydraulic forces are calculated using nonlinear hydraulic flow equations: 




h - ^ [-«0 - "2 " s <Pl - *2 ip] 


( 6 ) 
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2.3 (Continued) 

Piston forces calculated include the effects of viscous damping. Bp, and 
coulomb friction, F|:. For each actuator: 


The coefficient Cp is used to avoid a discontinuity at zero velocity. 

Force 


If linl ^ V > f**®" % ' 


IV 


If lipl < . then Cp = ^ 


'bw 






'bw 


( 8 ) 


Total hydraulic actuator forces and moments are then calculated for the 
equations of motion as follows: 





(9) 


*■ 


t- 
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2.3 (Continued) 



where the terms T. , T. , etc. are the terms in the transformation 
• ’ll ’l2 

from table coordinates to local actuator coordinates: 

h]=h7w (11) 

is the transformation from actuator coordinates to inertial 
coordinates. 

The servo electronics consist of actuator position and rate command 
signals and the electronic components shown in Figure 4. The use of the 
forward loop compensation network, the valve dynamics representation, and 
the position and rate feedback filter are optional. If B, for example, 
is input as a value less than unity, then the forward loop compensation 
network is not included in the simulation. Similarly, the valve dynamics 
and position feedback filter are neglected if < 1 and < 1, respec- 
tively. 

Table motion commands are input in the inertial coordinate system and are 
transformed to commands to the six actuators as follows: 

Define j^A 1 as the [A] matrix with the angles e, tj>, (}> replaced with the 
coirananded Euler angles e., i}>^. Then the commanded inertial components 

w V w 

i 
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Figure 4. Servo Electronics Block Diagram 
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2.3 (Continued) 
of actuator length are: 



and the correnanded inertial velocities of the table/servo attachment points 
are: 


where: 



(13) 






w. 



(14) 


Then the commanded actuator lengths and velocities are: 


"•c +r2 

i V "i ^i "i 


(15) 


1 









(16) 
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3.0 INPUT/OUTPUT DESCRIPTION 


3.1 INPUT DESCRIPTION AND PREPARATION 


Card 1 Format 13A6,A2 


TITLE 


TITLE - 80-character title to be printed on the top of the first 
c: page of output 

INERTIA AND GEOMETRY DATA 

Card 2 Format E12.6 


MT - Table mass 


Card 3 Format 6E12.6 
(INER(I), 1=1,6) 

INER - Table moments and procjucts of inertia 

(2) ' ^yy 

(3) - Izz 

(4) - 

(5) - . ■ 


Card 4 Format E12.6 


RXA - X table station of actuator swivel joints w.r.t. table c.g. 


■"4 S 
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3.1 (Continued) 

Card 5 Format 6E12.6 

((RfZA(J.I),0=l,2).I=l,6) 

RYZA(1,I) - y-table coordinate of I'th swivel joint w.r.t. 
table c.g, 

RYZA(2,I) - Z-table coordinate of I'th swivel joint w.r.t. 
table c.g. 

Card 6 Format 6E12.6 

{(YZF(J,I),J=1,2),I=1,6) 

YZF(1,I) - y-inertial coordinate of I'th floor swivel joint 
YZF(2,I) - Z-inertial coordinate of I'th floor swivel joint 

ACTUATOR STRUCTURAL DATA 


Card 7 Format 6E12.6 

ZETAE,MP.IAC,LC,LR,LO,EIR,LPM 

ZETAE " Actuator bending damping constant 


MP 

- Mass of rod and piston 


lAC 

- Moment of inertia of cylinder 

about floor swivel joint 

LC 

- Distance from floor swivel to 
seal at end of cylinder 

center line of piston rod 

LR 

- Length of piston rod 

• 

LO 

- Retracted length of actuator 


EIR 

- Bending modulus of piston rod 


LPM 

- Maximum stroke of actuator 

/ 


Note : The data above occupy two cards. 
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3.1 (Continued) 


FORCING FUNCTION DATA 


Card 8 Format 6E12.6 

(TCGC0(I),I=1,3),(TEA0(I),I=1,3) 

TCGCO - Initial inertial coordinates of table c.g. 

TEAO - Initial Euler angles of table coordinate system w.r.t. 
inertial system (rad.) 

Card 9 Format 414 

IM,NFREQ,IXF,NFFREQ 

IM - Mass matrix and geometry update option 
= 0, do not update after initialization 
= 1, update throughout time span 

NFREQ - Number of displacement frequency cases to run (max. = 18) 
(see note after Card 11) 

IXF - External force and moment option 
= 0, no external forces, or moments 

= 1, constant external forces and moments are to be applied 
to the table c.g. 

= 2, external forces and moments are sinusoidal 
NFFREQ- Number of external force and moment frequencies (max. = 18) 

Card 10 Format 6E12.6 
(0MEGAC(I),I=1,NFREQ) 

OMEGAC - Displacement command signal frequency (rad. /sec.) 

If OMEGAC(l) < 0, the command amplituc|es (Card 11) are 
assumed to be step velocities. 

Card 11 Format: 6E 12. 6 
(DLTAI(I),I=1,6) 
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3.1 (Continued) 


DLTAI - Sinusoidal amplitudes of X, Y and 2 commands for table c.g., 
displacement and table Euler angles, e, ^ and 4 > (if OMEGAC > 
0) or step velocities (if OMEGAC < 0) 

Note : Due to the interaction of Cards 10 and 11, it is logical that if 
one OMEGAC(I) is less than or equal to zero they must all be. 
Therefore, it is logical that NFREQ should be only one in that 
case. 


If IXF=0, skip the next two cards 

Card 12 Format 6E12.6 
(FMEXT(I),I=1,6) 

FMEXT - Magnitude of external forces and moments applied to the 
table c.g. 


If IXF=1, skip the next card 

Card 13 Format 6E12.6 
(0MEGAF(I),I=1,NFFREQ) 

0ME6AF = Frequencies of sinusoidal external forces and moments (rad. /sec.) 


HYDRAULICS DATA 


Card 14 Format 3E12.6 

PS,BETAE,KC 

PS ~ Supply pressure 

BETAE - Equivalent system bulk modulus 

KC - Valve pressure flovi coefficient 


Card 15 Format 6E12.6 
(CP(I),I=1,6),(BP(I),I=1,6) 
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3.1 (Continued) 

CP - Leakage coefficient across piston seals for each actuator 
BP - Actuator viscous damping coefficient for each actuator 
Note : The data above occupy two cards. 

Card 16 Format 4E12.6 

(A(I),I=1,2) (V0(I),I=1,2) 

A - Actuator push and pull stroke working areas 

VO - Initial hydraulic volumes of fully retracted actuator 

Card 17 Format 6E12.6 

(FF(I),I=1,6) 

FF - Coulomb friction force of each actuator 

ELECTRONICS DATA 


Card 18 Format 6E12.6 

(KG(I), 1=1,6) 

KG - Electronics and valve forward loop gain 

Card 19 Format 6E12.6 

(KF(I),I=1,6) 

KF - Displacement feedback and conmand gain 

Card 20 Format 6E12.6 

(KR(I),I=1,6) 

KR - Velocity feedback loop gain 

Card 21 Format 6E12.6 

(KPF(I), 1=1,6) 

KPF - Pressure feedback loop gain 
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3.1 (Continued) 

Card 22 Format 6E12.6 

(KRC(I), 1=1.6) 

KRC - Velocity command gain 

Card 23 Format 4E12.6 

ALPHA,BETA,(0ME6PF(I),I=1,2) 

I 

Break frequencies of | 

first order filters (rad. /sec.) 


KLrm - a 
BETA - e 

OMEGPF **’*,ci *^oco 


Card 24 Format 6E12.6 

ZETAS ,0MEGAS ,ZETAV ,0HEGAV 

ZETAS - Damping constant of second order filter on displacement 
and velocity feedbacks 

OMEGAS - Frequency of the displacement and velocity feedback filter 
(rad. /sec.) 

ZETAV - Damping constant of valve dynamics 
(HiEGAV - Frequency of the valve dynamics (rad. /sec.) 

INTEGRATION CONTROL DATA 

Card 25 Foraat 4E12.6,2I5 

TSTART ,TENU , (OUTFRQ ( I ) , 1=1 ,2 ) , IPROPl ,IPL0PT 
TSTART - Start time 

TEND - Stop time 

OUTFRQ(l) - Output frequency for printing (At, sec.) 

0UTFRQ(2) - Output frequency for plotting (at, sec.) 

21 
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3.1 (Continued) 

IPROPT - Print option • 

IPLOPT - Plot option 

Note : IPROPT and IPLOPT are of the form group 

number of the i*th group of data to be printed. These groups are 
explained in the next paragraph* Output Description. 

Cards 1 through 25 may be repeated as many times as desired. 

3.2 OUTPUT DESCRIPTION 

The output includes printed listings and plots of responses versus time. 
These responses are divided into four groups: 

Group 1 - Table Response Data 

a. Incremental inertial motions of the table c.g. 

b. Incremental angular motions 

c. Incremental velocities of the table c.g. 

d. Euler angle rates 

e. Table position errors 

Group 2 - Actuator Bending Data 

a. Bending frequencies of the actuators 

b. Y and Z lateral elastic displacements at cylinder rod seal 

Group 3 - Actuator Responses 

a. Actuator strokes 

b. Actuator velocities 

c. Actuator position error 

d. Net forces on the actuator pistons 
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3.2 (Continued) 

Group 4 - Complete Derivative. and Variable Arrays (may not be plotted) 


The variables and their derivatives (indexed by row) are listed below- 


Index Derivative Variable 



Variable Definition 


Table c.g. velocities 


Table rotational rates 


Lateral bending velocities of 
actuator along y. axis 

5 - 


Lateral bending velocities along 
z. axis 

94 


Inertial displacements of 
table c.g. 


Table Euler angles 


Bending deflections of actuators 
at top of cylinder 
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3.2 (Continued) 


Index 

31 

* 

ft 

ft 

36 

37 

ft 

ft 

ft 

42 

43 

ft 

ft 

ft 

48 

49 

I 

54 

55 

1 

60 


I * 

61 

1 


1 

66 

y- 


f i 

67 


1 


1 

f ■ 

72 


Derivative Variable 



Variable Def i ni t i on 


Bending deflections of actuators 
at top of cylinder 


"Push" hydraulic pressure on 
actuator pistons 


"Pull" hydraulic pressure on 
actuator pistons 


Derivatives of no-load valve 
flow 


No-load valve flow 


Filtered actuator feedback velocities 
from second order filter 


Filtei^ed feedback displacements 
from second order filter 


J 

\ 

I 


t 


I , 

: i* 


!l . 
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3.2 

(Continued) 


. 

Index Derivative 

Variable 

Variable Definition 

73 

• 



1 

1 1 

.1 



i 1 

1 

Voltage output of first order lag 

T 

i V 

• 

t 

filter in pressure feedback 

78 


a 



<’6 


- 

79 

i 



. - 

1 

, j- 

1 

Voltage output of high pass filter 

t 

f 

4 


84 

% 

% 1 

j in pressure feedback 

85 


1 

1 

1 

t 

i 

Voltage output of forward loop 

t 

T 

t i 

compensation filters 

90 


1 


3.3 

ERROR MESSAGES 




a. “THE INTEGRATION HAS FAILED AT T = XX.XXXX. ABORT AND 60 TO NEXT 
CASE." - occurs if for some reason the integration procedure cannot 
continue. The program will abort the case and attempt to process 
another one. 

b. "ACTUATOR XXX HAS STROKED OUT... ABORT AND 60 TO NEXT CASE." - occurs 
when an actuator has exceeded the maximum stroke. The program will 
abort the case and attempt to process another. 

c. "ERROR WHILE INVERTING MASS MATRIX... 60 TO NEXT CASE." - is self- 
explanatory. The user should look for errors in input which may 
cause a singular or ill-conditioned~mass matrix. 

d. There are also several error messages output by the integration routine. 
When one of these messages occurs, it is likely that an instability 

ha? occurred in the hydraulics or electronics caused by improper data. 
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4.0 OPERATING INFORMATION 

4.1 PROGRAM AND DATA SETUP 

The program may be input via standard EXEC II control cards on source or 
relocatable decks. All data for the program are input on cards. 

4.2 RUN INFORMATION ' 

Compilation time for the program is about 40 seconds. Representative runs 
of about one second of simulation time have averaged approximately 8 minutes. 
Run time depends greatly upon the frequency of the hydraulics and electronics 
inputs. 

5.0 PROGRAMING INFORMATION 

5.1 FLOW CHARTS 


Program flow charts are shown on the following pages. 
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ACAVAP: Actual Actuator Positions, Velocities, and Accelerations 
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ACDVC: Actuator Displacement, Velocity, and Acceleration Commands 

START 
AwOVC 


f START ^ 
V AwOVC j 


CALCULATE 

INERTIAL 

COMMANDS 


CALCULATE 
I' ,AH6. BODY 
, RATES AND 
ACCELERATIONS! 



CALCULATE 
tOMM. INERTIAL 
VELS. AND 
ACCELS. 


CALCULATE 
ACT. LENGTH 
INERTIAL 
COMPONENTS 


I ^ < 


I 


CALCULATE 
:OMM. ACTUATOR! 
.EN6THS, VELS, I 
AND ACCEL 
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ACIBD: Actuator Inertia and Bending Dynamics Parameters 
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DATARD: Data Input Routine 
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START 

DATARD 


READ & 
PRINT 
. TITLE 


im & PRINT ^ 
\ INERTIA & / 

Veometry / 

\ DATA / 


mu & PRINT ^ 
\ ACTUATOR y 
kTRUCTURAL/ 
\ DATA / 


im & PRINT, 
\ FORCING / 
\FUNCTION / 
\ DATA / 


JEAD & PRINT 
IrfYDRAULICs/ 
\ DATA / 


r ^ 


jEAD & PRINT, 

electronic/ 

\ DATA / 


OD & PRINT, 
k PROGRAM y 
\CONTROL 
\ DATA / 


RETURN 
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FCN: Derivative Evaluation Control Routine 




CALCULATE \ 
TRANS- \ 
FORMATIONS/ 
(TTRANS)/ 


I 


FORM 

MASS 

MATRIX 

(MMATRX) 


I 


INVERT 

MASS 

MATRIX 

(GJR) 




CALCULATE 
FORCES AND 
MOMENTS 
(CFM) 


I 


FORM 
COLUMN OF 
GEN. FORCES, 
(COLUM) 

“ 1 “ 
calculate 

FQNS. OF 
MOTION 
DERIVATIVES 
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Hydraulic and Servo Electronic Derivatives 











i 

I 
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INITIL: 
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Initializatlc.j Routine 
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OUTPUT: Print and Write Output on Drum for i-iots 






















{ 










D2-118544-2 


5.2 SYSTEM- STRUCTURE 

No overlo^y structure is required for this program. 

5.3 LIBRARY SUBROUTINES 

Several routines from the Johnson Space Center plotting library are used. 
These routines and their fu ictions are: 

FILMAV - advances the film 
6RDSET - sets line intensities 
GRID - forms a grid 
PLOTIV - plots an array of points 
PRINT - prints axis titles. 

These routines may not be available at another location; in which case, 
the best action would likely be to substitute routines with, the same 
ftinctions rather than obtaining these routines from the JSC library. 

5.4 PROGRAM LISTING 

A complete program listing is shown on the follov/ing pages. 
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5.5 SAMPLE PROBLEM • 

The input data and output, both printout and plots, for a sample case 
are presented on the following pages. In the example, the table is 
commanded to move sinusoidally in the y direction with an amplitude 
of 0.1 inch. Although 18 frequency cases are specified, output from 
the first frequency case only is presented. In this example the val‘*e 
dynamics, forward loop compensation filter, and position and rate feed 
back second-order filter are ignored. Rate command and rate feedback 
gains are set to zero. 
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